Supplementary Material: "Polarons in the radio-frequency 
spectrum of a quasi-two-dimensional Fermi gas" 
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In this supplementary material, we discuss the calculation of the radio-frequency spectra arising 
from confinement-induced dimers and polarons in a quasi-two-dimensional Fermi gas. We determine 
the dimer binding energies, including both the tight axial confinement and the nonzero transverse 
confinement. We provide the probabilities for dimer-to-dimer transitions and the shape of the dimer- 
to-scattering state spectrum. We also find the energy and quasi-particle weights for polarons in the 
two-dimensional gas and the corresponding spectra for polaron to polaron transitions. 

PACS numbers: 03.75.Ss 



We begin by reviewing briefly in §|T]the radio-frequency 
spectrum arising from confinement-induced pairs, in- 
cluding final state interactions, but ignoring many- 
body effects, using the method employed for the three- 
dimensional case by Chin and Julienne [l|. We consider 
mixtures of the three lowest hyperfine states of 6 Li, de- 
noted |1), 1 2), 1 3). For the conditions of our experi- 
ments in a 12 mixture at 720 G, the observed 2^3 
threshold spectrum is well described by a 12-dimer-to- 
13-scattering-state transition. In contrast, at 834 G, the 
predicted dimer spectrum is in marked disagreement with 
the data. In particular, we find that the difference be- 
tween the ground and excited state dimer energies is too 
small. In § |TT] we determine the energies for nonintcr- 
acting confinement-induced polarons. We find that the 
locations of the observed resonances for a 12 mixture near 
834 G are well modeled by the predicted energy difference 
between isolated state 2 polarons and state 3 polarons, 
in a bath of atoms in state 1. 



I. CONFINEMENT-INDUCED DIMERS 

A simple golden rule calculation gives the radio- 
frequency-induced transition rate out of the initial state 
to all possible final states Ri(uj r f) = ^2pRf<-i, where 
Rf^i = (2ir/h)\H' fi \ 2 5(E f - E t ~ hu rf ), with H' f . = 
USlfi(F\I) /2. Here, flfi is the Rabi frequency for chang- 
ing the hyperfine state of a single atom from the cho- 
sen populated state (i) to the initially unpopulated state 
(/) and (F\I) is the overlap between the initial and fi- 
nal wave-functions for the relative motion of the atom- 
pair. Since the center of mass energy does not change 
in the rf transition, Ef — Ei is the total change in the 
atomic hyperfine energy (= Tiw/j) plus the change in 
the energy of the relative motion of the pair Ep — Ej. 
Since £ F |(F|/)| 2 = 1, / dw rf Ri(u rf ) = (tt/2)^. We 
define a normalized spectrum I(ui) where Ri(ui r f) = 
(7r/2)ri 2 j I(uj) and w r / = Ufi + u, with u> the frequency 
relative to the (unshifted) free-atom hyperfine transition 
frequency. Then, I{lo) = £ F |(F|/)| 2 hS ( E F — Ei — Two) 



and j du)I{uS) = 1. 

To determine the spatial wavefunctions and the pair 
binding energies, we note that the range of the two-body 
interaction is small compared to the interparticle spac- 
ing as well as to the harmonic oscillator confinement 
scale l z = y/h/(mu z ). In this case, interactions be- 
tween atoms in two different spin states are well described 
by the s-wave pseudopotential in three dimensions 0], 
V(r) — (Airh 2 a/m) 6(r)d r (r...), where r is the distance 
between the atoms, m is the mass of a single atom and 
a is the magnetically tunable s-wave scattering length. 
The spatial wavefunctions are readily written in terms of 
the Green's function Ge{v) for the relative motion of the 
two atoms in the confining potential, which we take to be 
harmonic, with ground state energy Eq. The two-atom 
scattering states with energy E = E + E s , where E s > 0, 
take the form ip s (r) = ipp (r) — a Ge s (r) u' s (0), where 

ip^ is the input state and ip — u/r with u regular at 
r = 0. For the bound states, where there is no input, we 
have = — a GE b (j) u' b (0), where E = E - E b with 

Eb > 0. Using 9 r [r^f,(r)] r _>n = u' b (0) yields the equa- 
tion for the binding energy 0], 1 = — a d r [rGE b (r)] r _>0; 
where the right side projects out the regular part of G 
at r = 0. For a three dimensional harmonic trap, the 
Green's function is 



G e (r) 



dr,eiev n 



li) fM, (i) 



where E — Eg + ehcu z , with e = — eb for bound 
states and e > for scattering states. Here, fj{rj) — 

vWN"*)) with Pz = l, Px, y = Ux, y /u, and 
lj = h/\fWr After subtracting the irregular part of G 
and changing to real variables, we determine the dimer 
binding energies Ef, — ebhui z > from 



du 



V 47TW 3 



n t 



2/3,- u 
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(2) 



The dimer binding energy is significantly increased for 
nonzero transverse confinement. At resonance, where 



2 



l z /a —> 0, for v±/v z = we obtain Eb — 0.245 hv z , while 
for u±/v z = 1/25 we obtain E b = 0.290 hv z . At 842 
G in the shallowest trap, the binding energy of the 13 
dimer is increased from 0.15 kHz without transverse con- 
finement to 0.78 kHz with transverse confinement. We 
compute the pair binding energy Eb = eb TujJ z as a func- 
tion of magnetic field using the s-wave scattering lengths 
a measured in Ref. [3j. 

The scattering state is determined using u' s (0) = 

?4 0) (0)/{l + ad r [rG E .(r)] r -x)}. Using the relation 
(one-to-one correspondence) between the scattering 
length and the bound state energy, we have u' s (0) = 

u' { °\o)/{a [G Es (r) - G Eb (r)] r -x)}, where E b is the bind- 
ing energy corresponding to the scattering length a and 
we have used d r {r[G Es (r) ~ G E b (^)]r^-o} = [Ga(r) - 
GE b (r)]r->0) which is regular at r = 0. Then, the scat- 
tering state takes the form 



i> 8 



G E3 (r)uf\0) 
[G Es (r)-G Eb (r)} 7 



(3) 



where u'^\o) = i/jg(0), since the input state is regular 
atr = 0. 

For 2 — > 3 transitions in a 12 mixture at 834 G, the 
binding energies are small compared to the energy differ- 
ence between symmetric axial states 2hu> z , which are cou- 
pled by the s-wave scattering interaction. In this case, the 
Green's functions, and hence the 1-2 and 1-3 bound states 
and the 1-3 scattering states, are well approximated by 
the ground axial state component oc J dz<j>o(z) G tb (z, p), 
yielding the normalized bound state 



A b {z,p) = (j) {z)—j=K (Kp), 

\/7T 



(4) 



where <j>o(z) is the axial ground vibrational state for the 
relative motion, Eb — eb Hlo z is the pair binding energy 
for the given scattering length, and K 2 K 2 /m — Eb, i.e., 
k = ^fe~b/l z . The 1-2 bound state to 1-3 bound state tran- 
sition strength is then determined by the overlap between 
of the modified Bessel functions, -Ko( K i2/?) and -Kq(ki3p), 



12(13) 



The 1-2 bound to 



where Ki2(i3) is determined by e 
1-3 bound contribution to the lineshape is then given by 



hb{v) = e bb {q)S[is~ {e\ 2 



,13 



(5) 



where v z is the axial harmonic oscillator frequency in Hz 
and v is the rf frequency in Hz, relative to the bare 2- 
3 hyperfine transition frequency of ~ 83 MHz. For the 
axial ground state, the frequency integrated bound to 
bound transition strength can be written compactly in 
terms of q = ln(e£ 3 /e£ 2 ), 



dvlbbiy) = e bb {q) 



4sinli (g/2) 



(6) 



We plot eb b as a function of magnetic field in Fig. [T] for 
the general case, valid for both weak and tight binding of 
the 1 — 2 or 1 — 3 dimers, including the contribution of the 




600 700 800 

Magnetic field (G) 

FIG. 1: Bound 12 dimer to bound 13 dimer transition frac- 
tion versus magnetic field for v±/v z — 1/25. Upper curves 
(red online): Trap depth Uo = 280 /iK, v z = 82.5 kHz; Lower 
curves (blue online) (7o = 21 fjK and v z — 24.5 kHz. Solid 
lines denote the square of Frank-Condon overlap integral, in- 
cluding the first 50 even axial states. Dashed lines denote the 
corresponding results using only the ground axial state, Eq.(6] 



first 50 even axial states. For comparison, we show the 
weak binding approximation given by Eq. [SJ which in- 
cludes only the ground axial vibrational state. For small 
binding energies, increasing the trap depth significantly 
increases the pairing energy, increasing the overlap inte- 
gral for bound-to-bound transitions. 

From Fig. Q] we see that for a trap depth Uq = 21 pK 
at 720 G, the 2 — > 3 transition in a 12 mixture is dom- 
inated by bound-to-free transitions. In contrast, at 834 
G, the corresponding bound-to-bound transition is dom- 
inant and increases with increasing trap depth. 

In the same approximation, the two dimensional box- 
normalized (to area A) 1-3 scattering state takes the 
form ip e± (z,p) = <fro(z)ip e± (p), with E s = e ± huj z = 
H 2 k\/m. Assuming that the incident transverse state 
is the azimuthally-symmctric part (I = 0) of a box 
normalized plane-wave state exp(ik^ • xjJ/vM", i.e., 
Jo(k±p)/\/A, Eq. [3] gives 



AAp) = 



l 



Jo{k±p) 



ln(e & /ej_) + ni 



rH (k ± p) 



(7) 

We determine the overlap integral of the 1-2 bound state 
with the 1-3 scattering state and integrate the transition 
rate using density of transverse states, A/ (2Tr)2irk±dk±, 
to obtain the 1-2 bound to 1-3 scattering state contribu- 
tion to the lineshape, which takes the form of a threshold 
function, 
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(8) 
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Note that J dv I b f{y) = 1 — Cbb(q), as it should. We find 
that Eq. [5] well fits the radio frequency spectra obtained 
at 720 G, where the dimer binding energy is larger than 
the transverse Fermi energy. 



II. POLARONS IN A TWO-DIMENSIONAL 
FERMI GAS 

For radiofrequency 2 — » 3 spectra obtained near 834 G 
in a 12 mixture, we find that the difference between the 
calculated dimer binding energies significantly underes- 
timates the observed frequency shifts, as shown in Figs. 
2, 3, and 4 of the main paper. We consider the possibil- 
ity that the spectra may arise from transitions between 
polaronic states, as the polaron energy is significantly 
more negative than the corresponding dimer energy and 
is therefore energetically preferred 

To estimate the polaron energies, we consider either an 
isolated impurity atom in state 2 or in state 3, immersed 
in a bath of atoms in state 1. We employ the method 
described for polarons in three dimensions in the supple- 
mentary material of Schirotzek at al. @, which is based 
on the zero momentum polaron wavefunction proposed 
by Chevy Q, which for a polaron in state i — 2,3 takes 
the form 



|B ( )= ¥ J0i|0) i |FS) 1 + ^ Vkq |q-k) i 4 1 Cqil^i- ( 9 ) 

q<k F <k 

Here, the first term describes an impurity i of zero mo- 
mentum in a Fermi sea of atoms in state 1 for which the 
net momentum is zero. Collisions between the impurity 
and the background atoms couple the zero momentum 
impurity state to that with momentum q — k, producing 
a particle-hole pair from the Fermi sea of atoms in state 
1 with net momentum k — q, conserving the total zero 
momentum. 

For the 2D calculations, we replace the box normaliza- 
tion volume in the supplementary material of Schirotzek 
at al. [5| by the corresponding area A, so that the polaron 
energy in 2D takes the form 



q<k f 



where 



f- 1 (E pi ,q) = 



(Jo 



A 



k>k f 



Eq-k 



-E„ 



(10) 



(11) 



pi 



Here, tk = h 2 k 2 / (2m) and kp is the local Fermi wavevec- 
tor with Ep — Ti 2 kp/(2m) the corresponding local trans- 
verse Fermi energy. 

Following Zollner et al., Ref. we assume that the 
effective bare interaction U arises from a short range 
2D potential, so that the matrix elements Z7kk' = go /A 
are momentum independent, go can be rewritten using 



the physical two-body T-matrix element in 2D, 1/go — 
l/T 2B (k ) - (1/A) £ k l/(2e fc - 2e fco ). This method is 
similar to that employed previously d, 0, H|- Here, 
we choose the T-matrix element Tk'ko = T2bQ*-o)/A so 
that the scattering rate calculated using the generalized 
Golden rule reproduces the scattering rate obtained from 
the 2D flux corresponding to Eq. [7J Then T2s(ko) = 
h 2 f(ko)/m, with /(fc ) = 4Tt/[m + ln(E b /e±)]. Here, 
e± = 2e/c and Eb = ebh(jj z is the dimer binding energy 
calculated from Eq. [2J Both the Tri term and the fco de- 
pendence in T2B (ko) are canceled by corresponding terms 
in the sum (1/A) ]£ k l/(2e fc - 2e fco ), so that f^ 1 (E pi ,q) 
is independent of fco as it should be. For attractive po- 
larons with energy E p i < 0, we then obtain the simple 
integral equation e(Lj) = S(ij,ej), where 



-2du 



-L + ln[y/(l- f) 2 -u+(l 



I)] 



(12) 



Here, e(Li) = ti = E p i/Ep and Li = ln(El l / Ep) for an 
impurity in state i. Eq. [12] yields the polaron energies 
E p i = £j Ep for the initial and final states i = 2, 3, using 
the binding energies El 1 — el l huj z determined from Eq.[2] 
Fig. [5] shows the attractive polaron energies obtained 
from Eq. Q21 which agree with those obtained in Ref. [4j . 
We see that the polaron energy E p is a large fraction of 
the local Fermi energy, which can be much larger than 
the corresponding dimer binding energy Eb- 
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FIG. 2: Polaron energy E p /Ef versus \u[Eb/ Ef\, where Eb 
is the dimer binding energy and Ef is the local Fermi energy. 

For radiofrequency transitions between impurity states 
2 — > 3 in a bath of atoms in state 1, the momentum of 
the impurity does not change. We therefore assume the 
coherent part of the spectrum is given by 



I(fiuj) = Z 2 Z 3 S[huj - E p3 + E p2 ], 



(13) 



where | ifQ 3 ip 2 1 2 = Z 2 Z 3 is the square of the overlap 
integral between the part of the initial and final polaron 
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FIG. 3: Polaron quasiparticle weight Z versus 1u[Ei,/Ef], 
where Eb is the dimer binding energy and Ef is the local 
Fermi energy. 



states that yields the coherent part of the spectrum. We 
determine Z 2 = 1 9^02 1 2 using Z^ 1 = 1 — dT,(L 2 , e)/de, 
with e — » E p2 /E F , as described in Ref. @ for the 3D 
problem, and similarly for Z 3 = 1 9^03 1 2 - Fig. [3] shows 
the quasiparticle weight Z obtained from Eq. [12J For the 



dimer binding energies in a 12 mixture near 834 G, Table 
1 of the main paper, we find that Z 2 — 0.85 and Z 3 ~ 
0.94 for the shallowest trap depth, both close to unity. 
Hence, we expect that the overlap between the initial and 
final polaron states is strong and that transitions between 
polaron states should make an important contribution to 
the spectrum. 

In the limit that the dimer binding energy is small 
compared to the local Fermi energy, i.e., E b << Ep, 
one verifies from Eq. [T^] that the corresponding polaron 
energy yields the limiting form Q 



En 



-2E F 



]n(2E F /E b y 



(14) 



which can be interpreted as an effective mean field shift, 
since Ef is proportional to the 2D density n± (p) . At 842 
G, where the dimer binding energy is reasonably small 
compared to the transverse Fermi energy in our experi- 
ments, this formula overestimates the magnitude of the 
polaron energy difference E p ^—E p2 by about 10% for the 
most shallow trap. At 811 G, it overestimates the energy 
difference by about 40%. As E b is not small compared to 
E F for most of the data, it is not surprising that Ep can- 
not be adjusted in Eq.[l4]to give the measured frequency 
differences. 
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We measure radio-frequency spectra for a two-component mixture of a 6 Li atomic Fermi gas in 
the quasi-two-dimensional regime. Near the Feshbach resonance, where the transverse Fermi energy 
is large compared to the confinement-induced dimer binding energies for the initial and final states, 
we find that the observed resonances do not correspond to transitions between confinement-induced 
dimers. The spectrum appears to be well-described by transitions between noninteracting polaron 
states in two dimensions. 

PACS numbers: 03.75.Ss 



Quantum degenerate atomic Fermi gases, with mag- 
netically controlled interactions, are ideally suited for 
exploring pairing interactions in reduced dimensions [l|- 
[f|. Near a broad collisional (Feshbach) resonance, a 
two-component gas in three dimensions can be continu- 
ously tuned from a Bardeen-Cooper-Schrieffer (BCS) su- 
perfluid, which exhibits weakly bound Cooper pairs, to 
a resonant strongly interacting superfluid and finally to 
a Bose-Einstein condensate (BEC) of molecular dimers. 
Since 2002, degenerate strongly interacting Fermi gases 
have been studied in three-dimensional geometries, pro- 
viding a paradigm for strongly interacting systems in na- 
ture, from hig h temperature superconductors to nuclear 
matter [7H10|. In contrast to free space, where bound 
dimers exist only in the BEC regime, two-dimensional 
(2D) confinement also stabilizes bound dimers in the 
BCS regime P, HH . The interplay between confinement- 
induced pairing and many-body physics in 2D confined 
mesoscopic systems of several hundred atoms has not 
been previously explored and offers new challenges for 
predictions 0, d, Q . 

We study a quasi-two dimensional Fermi gas of 6 Li in 
the many-body regime. Radio-frequency (rf) spectra are 
obtained for a 50-50 mixture of the two lowest hyper- 
fine states denoted 1 and 2, by driving transitions to an 
initially empty hyperfine state 3 and measuring the de- 
pletion of state 2. At 720 G, well-below the Feshbach res- 
onance, the molecular 12 dimer binding energy is larger 
than the local Fermi energy, and the observed spectra 
exhibit the expected threshold form, arising from disso- 
ciation of 12 dimers into 13 scattering ("free") states. 
However, near the Feshbach resonance at 834 G the pre- 
dicted spectra for transitions between dimer states are 
in marked disagreement with the measurements, where 
the transverse Fermi energy is larger than the 2D dimer 
binding energy. Here, we find that the resonance lo- 
cations appear to be described by transitions between 
noninteracting polaron states, which describe an impu- 
rity atom in state 2 or in state 3, immersed in a bath of 
atoms in state 1. In this regime, polarons are expected to 



be energetically more favorable than the corresponding 
dimers 

Prior experiments have been performed with a dilute 
gas, where the local transverse Fermi energy /i is small 
compared to the binding energy Eb of a single dimer. In 
this case, the dimers are small compared to the interpar- 
ticle spacing and one expects the measured rf spectra to 
be consistent with predictions based on molecular bind- 
ing in the BEC regime and confinement-induced dimers 
in the BCS regime, as observed previously by measuring 
the threshold for free-to-bound transitions [5] and bound- 
to- free transitions [f|. 

Our experiments are performed instead in the many- 
body limit, where the the transverse chemical potential 
(i > Eb over most of the trap, but fj, ~ l.bhv z is small 
enough to be in the quasi- two-dimensional regime Q. 
For a simple BCS approximation in two dimensions, the 
trap-averaged rf transition rate to excite an atom from 
one populated state in a 50-50 mixture to an unpopu- 
lated noninteracting final state is cx L d 2 x.±\A\ 2 0[hu> + 
fi — '/i 2 + |A| 2 ]/w 2 , where u> is the radio frequency rel- 
ative to the bare atomic transition frequency and A is 
the pairing gap for the initial mixture. Here, we assume 
for simplicity that the temperature T « fi over most 
of the trap, so that the number of excitations is neg- 
ligible. For fj, » Eb, predictions for a 2D gas [l[ give 
A = \J1\x Eb in a mean field (BCS) approximation, yield- 
ing a threshold for the bound (Cooper pair)-to-free spec- 
trum at Tiuj ~ A 2 /(2/i) = Eb, which is just the dimer 
binding energy, as noted previously [gj]. Hence, the pre- 
dicted spectra are identical with the dilute gas limit, in 
contrast to our data taken in the same regime. 

Table Q] lists the experimental parameters for sev- 
eral different magnetic fields and trap depths. The 
CO2 laser standing-wave trapping potential is charac- 
terized by using parametric resonance in the weakly in- 
teracting regime near 300 G to determine the oscilla- 
tion frequencies of the atoms in the transverse direc- 
tions (y x ,Vy) and in the tightly confined axial direc- 
tion (v z ). The axial trapping potential is taken to be 



2 



B(G) 


Uo(fiK) 


i/ s (kHz) 


Nsite 




El 2 


El 3 




719 


27.5 


26.0 


1298 


1.87 


145 


2.9 


0.27 


809 


23.5 


24.0 


1951 


2.12 


12.6 


0.88 


0.53 


810 


294 


85 


1134 


5.73 


33.1 


6.69 


0.79 


832 


24.4 


24.5 


1620 


1.97 


7.25 


0.81 


0.66 


832 


274 


82.0 


1500 


6.36 


23.9 


5.79 


0.83 


832 


742 


135 


1800 


11.47 


39.1 


12.45 


0.89 


831 


1304 


179 


1250 


12.65 


51.9 


18.9 


0.91 


842 


24.4 


24.5 


1420 


1.85 


5.91 


0.78 


0.70 


842 


277 


82.5 


1617 


6.64 


21.4 


5.66 


0.85 



TABLE I: Parameters for the radiofrequency spectra: Trap 
depth Uo; Axial frequency v z ; Total number of atoms per pan- 
cake trap Nsite; Transverse Fermi energy Efj_ = hvxs/ Nsitu] 
El 2 and £ 6 13 are the dimer binding energies in kHz, for 
v±/Vz = 1/25. The dimer binding energies are obtained us- 
ing the Green's function method described in the supplemen- 
tary material [l5||. tbb is the corresponding bound dimer to 
bound dimer transition fraction (area under the bound to 
bound spectrum) determined from the binding energies ob- 
tained with the scattering lengths measured in Ref. [la |. 



tion fraction ebb = 0.27, we expect bound-to-free tran- 
sitions to dominate the resonance spectrum (TBJ. The 
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FIG. 1: RF spectra in a 1-2 mixture for a 2 — > 3 transitions 
in a quasi-two-dimensional 6 Li Fermi gas near 720 G. The left 
resonance occurs at the bare atomic transition frequency and 
the threshold resonance spectrum on the right is in very good 
agreement with predictions for molecular dimers. 



U axial — Uo sin 2 (2irz/ A) , where the trap depth is read- 
ily determined from v z in the harmonic approximation, 
U = m{is z \) 2 /2. We find u ± = ^PT/ = v z /25. The 
number of atoms per site is estimated using the lattice 
spacing of 5.3 /im and the number of atoms in the central 
part (along z) of the cloud, as measured by absorption 
imaging. 

To compare the measured spectra to predictions, we 
begin by determining the 2D dimer binding energies, 
E l b 3 = e l b 3 hv z for atoms in states i and j, as described 
in the supplementary material [H| and given in Table |U 
The calculation includes the finite transverse confinement 
of the trapping potential, where v_\_/v z — 1/25, which sig- 
nificantly increases the dimer binding energy, especially 
for weakly bound dimers. 

The contribution to the spectrum from a dimer-to- 
dimer transition is determined by computing the corre- 
sponding fraction e 00 . In the weak binding approxima- 
tion, including only the axial ground state part of the 
dimer wavefunction, we obtain, 



dvlbbiy) = e b b{q) 



4sinh^(g/2) 



(1) 



where q = ln(e£ 3 /e^ 2 ) for a 2 — >• 3 transition in a 12 
mixture. We plot ebb as a function of magnetic field in 
the supplementary material fl5j . There we show that 
similar results for ebb are obtained including the first 50 
even axial states, which determines ebb for both weak and 
tight binding of the 1 — 2 or 1 — 3 dimers. 

Fig.[T]shows the measured spectra at 720 G, well below 
the Feshbach resonance, where the molecular dimer bind- 
ing energy is larger than the local Fermi energy. At this 
field and trap depth, where the bound-to-bound transi- 



resonance near 150 kHz is well fit by a threshold func- 
tion for a quasi-two-dimensional gas using the calculated 
12 and 13 dimer binding energies (l5| . 
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where v is the rf frequency in Hz, relative to the bare 
atomic transition frequency, and E-, 
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e\^hv z is the 12 

dimer binding energy in Hz. Note that J dv Ibf(v) = 

1 — e bb{q), as it should. 

Next, we examine spectra for the 12 mixture near the 
Feshbach resonance, at 832 G, Fig. [2j at 810 G, Fig. El 
and at 842 G, Fig. |4j For several of the spectra, the 
bare-atom resonance peak on the left side exhibits a fast 
rise and a tail toward higher frequency, which we assume 
arises from a density-dependent mean field shift Jjj- At 
832 G, the sharp threshold remains at nearly the same 
frequency (which we have set to in Fig. [2]), even when 
the trap depth is increased to maximum, where the ideal 
gas Fermi energy at the trap center is Ep± = 12.7 /xK, 
i.e., ~ 260 kHz. We assume that the threshold loca- 
tion is determined by the lowest density region at the 
cloud edges, where the mean field shift is negligible. The 
observed bare-atom resonance frequency is in agreement 
with that obtained at high temperature as well as for 

2 — > 3 transitions with all the atoms initially in the 2 
state. 

Near 834 G, for the conditions of our experiments, the 
bound-to-bound transition fractions in Table U are close 
to unity [15j so that bound-to-bound transitions should 
dominate the spectrum above the bare atomic resonance 
frequency. To determine the frequency shift Ais between 
the bare atom resonance and the second resonance peak 
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FIG. 2: RF spectra for a 1-2 mixture (12 — > 13 transi- 
tion) near the Feshbach resonance at 834 G, versus trap 
depth Uo- The atom number remaining in state 2 is shown 
versus rf frequency: a) Uo = 21 /iK, v z — 24.5 kHz; b) 
U = 280 fiK, u z = 82.5 kHz; c) U = 742 ^K, v z = 135 
kHz; d) Uo = 1304 fjK, v z = 179 kHz. Vertical lines show the 
measured bare atomic resonance position (solid) and the pre- 
dicted frequencies for confinement-induced dimers: bound-to- 
bound transition resonance hv = El 2 - El 3 (dashed) and the 
threshold for the bound-to-free transition hv = El 2 (dotted). 
According to Table |IJ the bound-to-free transition (dotted 
line) should make a negligible contribution to the spectrum. 



to the right of it, we fit a threshold line shape to the 
atomic peak, which is convolved with a narrow Lorenz- 
tian. In this case, the peak position is shifted to the right 
of the bare atom threshold location, which is determined 
from the fit. A second Lorentzian, with a larger width, is 
fit to the shifted resonance peak. The results for Av are 
given in Table |TT] and do not agree with the predictions 
for dimer-to-dimer transitions, hAv^i 
shown as dashed lines in Figs. 

We consider the possibility that polarons and not 
dimers determine the primary spectral features, so that 
the difference between the initial and final state polaron 
energies determines the observed frequency shifts. We 
assume that the coherent part of the spectrum is ap- 
proximately given by Z S[hu) — E p (3, 1) +E p (2, 1)], where 
E p (i,l) is the polaron energy for an impurity atom in 
state i = 2, 3 immersed in a bath of state 1. Since the mo- 
mentum does not change in the rf transition Z ~ Zi Z% is 
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FIG. 3: RF spectra for 2 — > 3 transitions in a quasi-two- 
dimensional 6 Li Fermi gas at 810 G: a) Uo = 23.5 /iK, v z — 
24.0 kHz; b) Uo = 294 fjK, v z = 85.0 kHz. Vertical lines show 
the measured bare atomic resonance position (solid) and the 
predicted frequencies for confinement-induced dimers: bound- 
to-bound transition resonance hv — Eb\2 — E\,\z (dashed) and 
the threshold for the bound-to-free transition hv — E^n (dot- 
ted). 



determined by the overlap of the initial and final polaron 
momentum space wavefunctions [lij . For experiments 
at 832 G, we find that the individual polaron quasipar- 
ticle weights are Z 2 > 0.8 and Z 3 > 0.9 15], assuming 
that the local density does not change in the transition. 
Therefore, we expect strong overlap between the initial 
and final polaron states, so that transitions between po- 
laron states should make an important contribution to 
the spectrum. 

We determine the 2D polaron energies for an isolated 
impurity atom in state 3 or state 2 in a bath of atoms 
in state 1, using the method described in the supple- 
mentary material [lij]. The method is based on that de- 
scribed for a 3D gas in the supplementary material of 



Schirotzek at al. 12|, which utilizes the zero momentum 
polaron wavefunction proposed by Chevy [3]. We ex- 
tend that method to the 2D problem , by renormalizing 
the interaction strength as in Refs. P,[H,[ll. Using the 
calculated dimer binding energies to determine the po- 
laron energies, we find that for an atom in state 3 at 832 
G, the energy for an isolated polaron is attractive and 
more than half of the Fermi energy, much larger than the 
corresponding dimer binding energy. For an impurity in 
state 2, for which the 12 scattering length diverges, the 
polaron energy is attractive and somewhat larger than 
the Fermi energy. In this case, the polaron is localized to 
approximately the interparticle spacing, but is not small 
compared to the interparticle spacing, as it is at 720 G. 



4 



l — 1 — i — 1 — i — 1 — i — ■ — i — ■ — i — 1 — i — 1 — i — ■ — i — 1 — i — 1 — i — ■ — i — ■ — i — 1 — i — 1 — i — ■ — r 




-10 10 20 30 40 50 60 

Frequency (kHz) 

FIG. 4: RF spectra for 2 — > 3 transitions in a quasi-two- 
dimensional 6 Li Fermi gas near 842 G: a) Uo = 24.4 fiK, v z = 
25.0 kHz; b) U = 277 fj,K, v z = 82.5 kHz. Vertical lines show 
the measured bare atomic resonance position (solid) and the 
predicted frequencies for confinement-induced dimers: bound- 
to-bound transition resonance hv — Eti2 — i?bi3 (dashed) and 
the threshold for the bound-to-free transition hv — Eti2 (dot- 
ted). 



In calculating the polaron frequency shift, we assume 
the peak position is determined by a local Fermi energy 
Ep = Ai Ep±, where Ep± is the ideal gas global Fermi 
energy given in Table |U Fixing Ai = 0.67, we obtain the 
polaron frequency shifts hAv po i aron = E p (3, 1) — E p (2, 1), 
which are compared to the measurements in Table [II] 
For the data just below and just above the resonance, 
at 810 and 842 G, Table ITT1 we find excellent agreement. 
Further, at 832 G, the calculated polaron frequency shifts 
agree very well with the measured frequency shifts at 
all four trap depths, Table UH demonstrating nontrivial 
scaling with the trap depth and axial frequency. 

We have also obtained spectra for a 50-50 mixture of 
states 1 and 3 near the Feshbach resonance at 690 G, 
where we drive either the 13 — > 12 transition or the 
13 — » 23 transition. In both cases, the binding energy of 
the final state 12 or 23 dimers is large compared to the 
local Fermi energy and suppresses the bound-to-bound 
transition probability. For the high density regime of the 
experiments at the lowest temperatures, the bound-to- 
free part of the spectrum is not separated from the bare 
atom transition peak, and probably contain a dimer to 
free contribution, while at higher temperatures, there ap- 
pears to be a polaron peak. These spectra are consistent 
with the predictions for polaronic transitions. However, 
more theoretical work is needed to explain the detailed 
shapes of all of the measured spectra, which may be im- 
proved in future experiments by local measurements @. 
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TABLE II: Frequency shift Av between the bare atom peak 
and the second resonance peak for a 12 mixture near the Fes- 
hbach resonance at different trap depths. The corresponding 
axial trap frequency is v z . The measured values of Av are 
compared to the values calculated assuming a transition from 
a polaron in state 2 to a polaron in state 3, in a bath of atoms 
in state 1. For determination of the polaron frequency shift, 
the dimer binding energies are taken from Table U and the 
local Fermi energy is taken to be Ef = 0.67 Ep±, where Ep± 
is the ideal gas global Fermi energy given in Table U 



In summary, we have observed rf spectra in a dense 
quasi-2D Fermi gas regime near a Feshbach resonance, 
where the dimer binding energies for both the initial and 
final states are smaller than the transverse Fermi en- 
ergy. The spectra appear to be explained by transitions 
between noninteracting polaronic states, despite the 50- 
50 mixture employed in the experiments. These results 
support the conjecture [Hj that at low temperature, a 
strongly interacting Fermi gas in a balanced mixture of 
two spin states may be described as a gas of polaron 
pairs. 
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